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The Radau-Darwin Approximation in the Theory of the 
Figure of the Earth. 

By K. E. Bullen, University of Melbourne. 


[Received by the Editor , September 2, 1943; issued separately , March, 1944 .1 


In a recent estimation (1942) * of the ellipticities of surfaces of 
equal density within the Earth, the writer made use of an approxi¬ 
mate equation first used by Radau (1885) and Darwin (1900) in 
problems on the theory of the figure of the Earth. Some further 
argument is necessary to justify the use of this approximation. 

Let c be the ellipticity of the surface of equal density (within 
the Earth) whose mean radius is r, m the mass of the matter en¬ 
closed by this surface, and p the density at points of this surface. 
For present purposes it is sufficiently accurate to assume that the 
distribution of mass within the Earth is consistent with a hydro¬ 
static state. The equation 


d 2 c 8?rpr 2 dc j 1 _ ^pr 3 | 6^ _ 

dr 2 + m dr ( m / r‘ J 

is then found (see Jeffreys, 1929) to give the variation of e with r 
to sufficient accuracy. 


When the density variation is known to a sufficient degree of' 
precision, as is the case now (see Bullen, 1940 and 1942A), it is 
theoretically possible by numerical integration of (1) to determine 
c as a function of r throughout the Earth. Direct use of the equation 
(1) would, however, involve much heavy labour, and this can be 
avoided by the use of the Radau-Darwin approximation, to be de¬ 
scribed shortly. But a point of logic arises in connection with the 
use of this approximation in that the approximation itself is not 
valid unless the true values of e obey certain restrictions. Actually 
the values of c as found in Paper I are such as to obey these re¬ 
strictions; but since these values were calculated using the approxi¬ 
mation, it would be arguing in a circle to infer from this result alone 
the validity of the approximation. The purpose of the present paper 
is to give the additional argument necessary to show that the Radau- 
Darwin approximation, and hence also the various results that its use 
entails, are indeed valid. The "writer has considered this point previ¬ 
ously (1936), but the following argument is more precise. 

As a preliminary to making his approximation, Radau modified 
the equation (1) by introducing parameters y and p c , where 

r dc 


y = 


dr 


( 2 ) 


and p 0 is the mean density of tli§ matter within the equal-density 
surface of ellipticity e; thus, to a sufficient approximation, 


* This paper will be hereinafter referred to as Paper I. 
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Transactions . 


Po r3 = J* * Spi^dr, (3) 

and m = ^ Po r 3 . (4) [ 

Replacement of € and m in (1) by 77 and p 0 leads (without further | 

approximation) to “ Radau's equation ” 

^|por 5 (1 + */)*[ = 5 Po r 4 if ( 1 /), (5) 

where ^ ( 77 ) = (1 + £ 77 — T h v 2 ) (1 + 77 ) ( 6 ) 

Use of the equation (5) would give sufficiently accurate values 
of e, but the process of calculation would still be severe. Radau's 
approximation was to put $ ( 77 ) equal to unity in (5) ; and rests on 
the circumstance that ^ ( 77 ) is equal to unity within 8 parts in 
10,000, if 77 lie& in the range 0 ^ 77 ^ 0*6. If, therefore, it can be 
established that 77 does in actual fact lie within this range inside 
the Earth, the approximation is justified. 

It will be convenient to re-write (5) in the form 

f o P ° r4 $ Of) dr = i po 1 * 5 (1 + 7 /) i * ( 7 ) 

The corresponding approximate equation is then 

J*' p 0 r 4 dr = i Po r 3 (1 + v )i, ( 8 ) 

the latter being equivalent to the equation (3) of Paper I, which 
was the basis of the writer's recent calculation of values of e. This 
equivalence may be readily inferred from Darwin's work on 
the theory of the figure of the Earth, or deduced in the present 
notation as follows. As in Paper I, introduce z where 


zmr 2 = I, (9) 

I being the moment of inertia of the matter inside the equal-density 
surface of ellipticity c. Then since, approximately, 

1 — J l f irpr 4 dr, ( 10 ) 

we have by (4) and (9) 

z Po v B = 2 J** pr 4 dr. (11) 

From (3) we can deduce that 

3 J"*pr 4 dr = p»r 5 — 2 p 0 (r)r 4 dr. ( 12 ) 

From (11) and ( 12 ), .we obtain 

(1 — V ) Por 8 = 2 f T o Po r 4 dr. (13) 

Finally, by ( 8 ) and (13) we obtain 

1 _ V = » (1 + V)K ' (14) 


which agrees with equation (3) of Paper I. Hence use of this 
equation (14), as was made in Paper I, is equivalent to using ( 8 ) 
to calculate values of 77 and c. 

The validity of the approximation in question will therefore be 
established if we can show that the use of the approximate equation 
( 8 ) in place of (7) introduces sufficiently small errors in 77 . We 
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shall now use rji to denote values of 77 as determined using ( 8 ), and 
shall keep 77 to denote the more accurate values that satisfy (7). 
Thus we re-write ( 8 ) as 

f: p»r 4 dr = i p 0 r 3 (1 + Vl )K (15) 

From the fact that the density increases steadily from the 
Earth's outer surface to the centre it can be shown using (1)—see 
Jeffreys (1929)—that e must increase steadily from the centre to the 
outer surface. Hence by ( 2 ), rj cannot be negative at any point 
within the whole Earth. Now the maximum value of ^ ( 77 ), as 
given by ( 6 ), for 77 ^ 0, is found to be 100074. Hence by ( 7 ) 

i PoY 5 (1 + V)* < 1-00074 fl p 0 r 4 dr. 

Hence by (15), 

(1 + 17 )* < 1 00074 (1 + Vl )*; 
i.e., 1 + V < 10015 (1 + 77 X ); 

i.e., V — Vi < 0 0015 (1 + Vi), (16) 

It follows that if the values of 77 as calculated using ( 8 ) need a 
positive correction, this correction will be less than 0.0015 (1 + 771 x 
for all values of r. Now the maximum value of 771 as determined in 
Paper I is 0*56. Hence the maximum positive correction to be 
applied to 771 to give 77 is less than 0*003. Thus if any of the previ¬ 
ously calculated values needs a positive correction, this will be such 
that 77 still lies in the range 0 ^ 77 ^ 0 * 6 . 

If on the other hand, '771 needs a negative correction for any 
value of r, the corresponding value of 77 will still lie in this same 
range, since as already pointed out, 77 cannot be negative. 

As we have thus shown that throughout the Earth 77 must lie 
in the range 0 ^ 77 ^ 0*6, it follows that the Radau-Darwin 
approximation is valid. This conclusion, moreover, substantiates the 
reliability of the calculations in Paper I, as well as the various other 
applications of the approximation, including notably Darwin’s esti¬ 
mation of the connection between the moment of inertia of the 
Earth and the ellipticity of its outer surface. 
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